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L'outil statistique en biologie du sol 
IX. Analyse spectrale et distribution spatio-temporelle 
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Synopsis: The technique of spectral analysis, described in some detail, and a new 
smoothing method were applied to the analysis of the spatio-temporal distribution of soil 
oribatid mite populations. 
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INTRODUCTION 


La variation des effectifs dans les recensements faunistiques du sol est à la 
fois temporelle et spatiale. D'habitude, les écologistes s'intéressent à exploiter plus 
en détail l'information temporelle tandis qu'ils s'efforcent de minimiser la variation 
spatiale (hétérogénéité spatiale) par l'intermédiaire des techniques d'échantillonnage 
(CANCELA DA Fonseca & VANNIER 1969), où ils sont moins intéressés par les 
variations autour des moyennes que par les moyennes elles-mêmes (PLATT & 
DrNMAN, 1975). Or, les écologistes tendent à se servir séparément de ces deux 
composantes de la variation totale comme s'il s'agissait de phénoménes indépen- 
dants; il en résulte des pertes d'information (PLATT & DENMAN, 1975). 


En ce qui concerne les outils analytiques en écologie et biologie du sol, la 
technique des corrélogrammes a été quelque fois utilisee dans l'analyse de la 
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structure de séries chronologiques (LEBRUN, 1971; BRUNEL & CANCELA DA FONSECA, 
1977; BELLIDO, 1979; DELETTRE & CANCELA DA FONSECA, 1979; CANCELA DA FONSECA, 
1980) et la technique des indices d'hétérogénéité ou d'agrégation utilisée assez 
souvent quand il s'agissait des questions de distribution spatiale (DEBAUCHE, 1958, 
1962; Воммет, 1964; BERTHET & GERARD, 1965, 1970; LEBRUN, 1965; CANCELA DA 
Fonseca, 1966; Lioyp, 1967; Usher, 1976; CANCELA DA Fonseca & SrAMOU, 1982; 
entre autres). 


Par contre, la méthode de l'analyse spectrale n'a pas été beaucoup utilisée, 
non seulement en écologie du sol, mais également dans les autres branches de 
l'écologie (voir Hacker et al, 1973; PLATT & DENMAN, 1975; SHUGART, 1978; 
LEGENDRE & LEGENDRE, 1979; Orp, 1979; RirLEY, 1981), mis à part l'écologie 
marine (par exemple : PLarr et al., 1970; IBANEZ, 1972; PLATT, 1972; LAUREC & 
BLANC, 1974; Denman & PLATT, 1975; EcoNoMiDis et al., 1988) et l'écologie végétale 
(dont HiLL, 1973; Оѕнев, 1975; RirLEv, 1978; KENKEL, 1988). Bien que l'analyse 
spectrale ait été utilisée principalement dans l'étude des séries chronologiques 
(BLACKMAN & Tukey, 1958; Hannan, 1960; Cox & MILLER, 1965; JENKINS & 
Warrs, 1968; DAVENPORT, 1970; ANDERSON, 1971; CHATFIELD, 1975, 1980; Вох & 
JENKINS, 1976; KENDALL, 1976; BARTLETT, 1978; AZENCOTT & DACUNHA-CASTELLE, 
1984), elle а pu être aussi employée dans l'analyse des distributions spatiales 
(BARTLETT, 1975, 1978; BENNETT, 1979; RipLEv, 1981). 


L'analyse spectrale, qui est un cas particulier d'analyse de variance d'une série 
de données (PLArr & Denman, 1975), données soit temporelles, soit spatiales. 
constitue cependant un outil important destiné à rendre apparentes les caractéristi- 
ques latentes de la structure de cette série de données, en particulier les périodicités 
difficilement décelables à l’œil nu (Jacos, 1980). Par cette analyse, la variance de 
la série autour de sa moyenne arithmétique est séparée en des contributions à des 
fréquences qui sont des harmoniques de la longueur de la série (PLATT & DENMAN, 
1975). 


Dans ce travail, nous nous proposons de décrire les techniques spectrales et 
de les appliquer par la suite, à titre expérimental, à l'analyse de l'information spatio- 
temporelle apportée par des recensements des populations d'Acariens Oribates 
édaphiques. La cinétique de populations recensées périodiquement au long de 
transects, ordonnées donc dans le temps et dans l'espace, est comparée à celle de 
populations également recensées périodiquement, mais au hasard. Notre hypothése 
est : Peut-on déceler, outre la variabilité temporelle, une variabilité spatiale différen- 
tielle au long des transects? 


1. — MATÉRIEL 


Les techniques spectrales ont été appliquées à l'analyse de la variation spatio-temporelle des 
populations de l'Acarien Oribate Achipteria holomonensis Cancela et Stamou (CANCELA DA FONSECA & 
SrAMOU, 1987). Les données analysées ont été obtenues par le recensement annuel de ces populations 
sur 12 transects de 1 m de largeur séparés par des lignes parallèles, recensement effectué d'octobre 1978 
a novembre 1979 dans une chênaie de la forêt du Mont Holomon, Halkidi, Grèce (STAMOU, 1981). 
L'Oribate était extrait, par la technique de Berlese-Tullgren, de la litière prélevée autour des rhizosphéres 
des graminées Festuca sp. (F) et Arrhenaterum elatius Mert. et K. (A) qui dominaient dans la parcelle 
expérimentale (50 x 15 m). Chaque fois, dans un des 12 transects, pris au hasard, cinq prélévements 
contigus circulaires de Ø 5 ст (cf. CANCELA DA FONSECA et al., 1967) étaient effectués autour de 
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chacune des six rhizosphéres, équidistantes et ordonnées sur un transect linéaire: F—-A—F—A—F—A 
( fig. 1). soit un total de 30 prélèvements par transect (tab. 1) 
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FIG |- Schema d'echantillonnage d'Achipreria hülomanensis autour des rhizosphéres des graminées елиса sp. (F) et 


Arrhenatherum clarius (AY 


Ces mêmes techniques furent appliquées à l'anaiyse de l'information spatio-temporelle d'un recense- 
ment annuel d' Achipteria coleoptrata (Linné) effectué de mai 1971 à avril 1972 dans les couches L et F 
de la litière de la hêtraie de la réserve biologique de la Tillaie, forêt de Fontainebleau ( ATHIAS-HENRIOT 
& CANCELA DA FONSECA. 1976). Chaque mois, huit prélèvements de litière, de mêmes dimensions que 
les précédents, étaient pris au hasard dans la parcelle expérimentale (tab. 11). Les prélèvements étaient 
numérotés de | à 8 par rapport à leur distance à un des côtés, toujours le méme, du parallélogramme 
expérimental. 


Les données furent disposées en deux types de tableaux : 


(a) Nombre d'animaux par ргёіеуетепі : Dates x Rhizosphères x Prélévements (tab. 1. А) ou 
Dates x Prélèvements (tab. I) 
b) Nombre moyen d'animaux prélevé par rhizosphére - Dates x Rhizosphères (tab. 1. B) 


Afin d'analyser l'information spatio-temporelle des données contenues dans ces deux tableaux par 
l'analyse spectrale. les données ont été réarrangées de façon à être placées sur une seule ligne. la méme 
relation d'ordre ayant été choisie pour tous les transects. Par exemple : 


(Moss 1) (Mois 2) (Mois 3... 11) (Mois 12) 
FI-AI-F2-A2-E3-A3- — FI-AI-FE2-A2-F3-A3-  FI-Al-...-A3- FI-AI-F2-A2-F3-A3 


La question à laquelle il faut repondre est donc : « En groupes de combien de prélévements ou de 
rhizosphéres. la décomposition de la série de données apporte-t'elle une contribution importante à la 
variabilité totale de l'information spatio-temporelle? » 


П. — DÉFINITIONS ET CONCEPTS 
A) Analyse spectrale 


Soit une série discrète d'observations, X,(x,: 1=1,2, ... №). La série X, peut être 
considérée comme unc réalisation du processus aléatoire stationnaire X (1). La stationnarité 
du processus aléatoire X (1) implique que les distributions gaussiennes de probabilité associées 
restent invariantes par translation dans le temps, c'est-à-dire que les éléments de seconde 
ordre qui décrivent entiérement chaque processus stationnaire — sa moyenne, sa variance, 
sa fonction d'autocovariance et la transformée de Fourier de celle-ci — restent invariantes 
dans le temps (HuvsERECHTS, 1975). Nous supposerons que les processus à analyser sont 
toujours stationnaires. Or les fonctions d'autocovariance ou d'autocorrélation et leurs trans- 
formées de Fourier — les fonctions de densité spectrale — sont des techniques équivalentes 
pour décrire des processus stationnaires (CHATFIELD, 1975). 
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elatius (A) au long des transects F 1-A 1-F 2-A 2-F 3-A 3. Mont Holomon, Grèce. 


Transect 


Tas. I. — Effectifs d'Achiptería holomonensis dans les rhizosphéres de Festuca sp. (F) et d'Arrhenaterum 


A. Valeurs recensées 
|Moyenne 
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|Moyenne 
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Tas. I. — (suite). 


——— + ———— ————— 
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Tas. 1 (suite et fin). 


B. Valeurs moyennes 


Adultes 
Transect Mois Total Moyenne ÉcartT | 
| 1 2 3 1 5 6 7 8 9 10 п 12 
FI mv 26 24 0.8 9.0 5.0 3.2 3.8 44 0.8 18 Fa 2.0 37.0 3.08 2.30 
Al 6.0 1.0 1.6 1.6 8.4 2.0 3.8 1.8 1.4 1.0 2.6 8.0 39.2 3.27 2.70 
E aaa rx à 54 6.0 2.0 3.6 0.8 9.2 6.8 1.0 1.8 1.4 22 38 | 44.0 3.67 2.66 
Э. «венаи d 1.6 0.2 46 6.2 1.2 3.2 0.8 1.4 1.6 2.4 1.2 29.4 2.45 1.89 
ЕУ EERE RTE 6.0 4.0 L6 6.6 5.4 0.8 4.6 3.8 1.8 2.6 2.8 6.0 46.0 3.83 1.93 
АЗ. „на Бумага 5.0 3.4 1.0 5.6 4.8 44 1.0 0.8 1.6 1.4 1.4 5.6 360 300 1.97 
TOTAL, sas 300 184 э 3L0 30.6 20.8 232 126 8.8 98 126 266 
Moyenne . e| $00 307 120 5.17 5.10 3.47 386 2.10 147 163 210 443 
ШЕЕ, oes а 126 181 066 255 2.48 3.11 189 160 0.37 0.54 0.65 2.58 
Immatures 
Transect Mois Total Moyenne Écart-T | 
1 2 3 4 5 6 7 8 9 10 п 12 
2.6 1.8 24 180 21.0 0.4 1.4 4.6 1.6 0.8 1.6 0.6 56.8 473 7.01 
7.2 42 32 7.4 7.4 1.6 2.6 1.6 1.4 1.4 L8 10,8 50.6 4.22 3.18 
6.0 74 3.6 6.4 11.4 4.0 1.6 1.6 1.8 0.8 2.0 1.6 48.2 4.02 3.20 
6.2 3.0 0 44 25.6 1.0 0.8 2.0 3.0 0.4 2.2 1.2 49.8 4.15 17.22 
46 206 6.0 72 44 8.0 1.6 2.6 L8 0.6 6.4 2.6 66.4 5.53 5.31 
| 60 54 2.2 6.2 11.6 2.0 0.8 0.4 3.6 2.2 1.6 42 46.2 3.85 3:13 
316 424 174 496 81.4 17.0 88 128 132 62 156 210 
Moyenne. ..... | 5.43 7.07 290 8.27 13.6 2.83 1.47 2.13 220 1.03 260 3.50 
Écart-T. : 1.62 691 1.97 489 8.13 2.81 0.67 1.41 089 0.66 1.88 3.79 
TorAUX 
Transect Mois Total Moyenne Écart-T 
1 2 3 4 5 6 7 8 9 10 il 12 
5.2 42 32 270 26.0 3.6 52 9.0 24 2.6 2.8 2.6 93.8 7.82 8.92 
13.2 5.2 4.8 9.0 15.8 3.6 6.4 3.4 3.0 20 44 188 89.6 7.47 5.54 
114 134 5.6 100 12.2 13.2 8.4 2.6 3.6 22 4.2 5.4 92.2 7.68 4.24 
11.2 4.6 0.2 90 31.8 22 4.0 28 44 2.0 4.6 24 79.2 6.60 8.50 
10.5 246 76 13.8 9.8 8.8 6.2 6.4 3.6 32 9.2 8.6 112.4 9.37 5.62 
11.0 8.8 32 118 16.4 6.4 1.8 1.2 52 3.6 30 9.8 82.2 6.85 4.71 
62.6 608 246 806 1120 378 320 244 222 156 282 476 
Moyenne. ..... 104 10.1 4.10 13.4 18.7 6.30 5.33 423 370 260 4.70 793 
Ber T. lesus | 2.72 7.90 2.52 690 8.49 413 226 290 099 067 233 612 
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Les fonctions d'autocovariance et d'autocorrélation du processus stationnaire X (1) sont 
définies par 


N-h 


y UO) S OUN) Y. (ХХ (е) иј (h=0, 1, 2, ... N—1) (1) 


1=1 


p(A) = y (h)/y (0) y (0/02 (2) 


où h est un entier représentant le décalage ou retard, et u( = x) et с2(= 52) sont respectivement 
la moyenne et la variance totale des valeurs de la série de données, et où les valeurs de y (A) 
ne dépendent que du décalage h (Сох & Мишев, 1965; CHaTFIELD, 1975). 


L'analyse spectrale étant l'étude du processus aléatoire stationnaire X(r) du point 
de vue de son contenu fréquentiel (Jenkins, 1961), leurs fonctions d'autocovariance et 
d'autocorrélation admettent une représentation spectrale dans le domaine de fréquence (!), 
représentation traduite par les transformées inverses de Fourier respectives (Cox & MILLER, 
1965: CuarriELD, 1975, 1980), dans leurs formes exponentielles et trigonométriques : 


ve | вано [ e' f (o) do (3) 


J-5 -r 


-| cos wh frin | cos wh f (о) do 
et 


s»- | “кчө | e'*^ f* (w) do (4) 


-a 


-| cos share | cos wh f * (o) do 


ou 


Е* (w) =F (0)/02 


et où o est la fréquence angulaire, F (œ) la fonction de distribution spectrale ou spectre de 
puissance et F*(w) sa forme normalisée (CuarrigLD, 1975, 1980). Cox et Mitten (1965) 
définissent la forme normalisée de CuarrigLD (1975, 1980) comme la réelle fonction de 
distribution spectrale. Dans un cas comme dans l'autre, ces auteurs prennent pour base le 
théoréme de Wiener-Khintchine, les premiers considérant que l'existence de la fonction 
d'autocorrélation, et le troisième que l'existence de la fonction d'autocovariance dans les 
processus stochastiques stationnaires impliquent qu'il existe une fonction monotonique 
croissante qui est respectivement F* (w) ou F (o). 


(*) La fréquence. nombre de cycles par unité de temps, est de la forme : f — 0/(2 x); l'inverse de la fréquence 
=Т= 1//= 2 n/a — est la période ou longueur d'onde, & étant la fréquence angulaire en radiants par unité de temps 
(Камоли. 1976: CHATFIELD. 1980). 
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Tas. П. = Effectifs d' Achipteria coleoptrata dans huit prélèvements pris au hasard dans la litière 
de hétre (couches L et F). La Tillaie, France. 
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La dérivée de la fonction de distribution spectrale par rapport à la fréquence angulaire 


définie la fonction de densité spectrale ou spectre, f(«) ou sa forme normalisée f * (w) : 


ou 


f (e) = dF (w)/do 


f[* (o)—d F*(o)/do 


(5) 


(6) 


qui sont respectivement les transformées de Fourier de l'autocovariance et de l'autocorréla- 
tion [cf. équations (3) et (4) du processus aléatoire stationnaire X (1)] : 


F= Y, he =N У yth) cos oh 
gey һ=- ж 


(7) 
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Тав. П. = (suite et fin). 


Preievemenis 


1 
s w 
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E 
5 4 
^ 
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Tots 16 
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et 


/*(®)=[!/(2л)]) Y, pe "^-[1/27)] 


- ханы 


COUCHE F 
Immatures 
Mois 
4 5 6 1 
3 x 5 ау 
1 
2 i 
D 4 4 
х 
1 3 2 
1 4 2 ^ У 
z 6 4 
п R? и 3 ж 
213 050 135 275 475 
23 232 2941 212 ай 
Adultes 
Mois 
4 5 6 1 
i 
i 
1 1 ' 
П 
' П 
2 i 
1 
2 E E a 
025 025 uso oso 
045 оле 054 
Totaux 
Mais 
4 E [ 7 
3 * «ау 
1 1 
1 3 І ' 
1 a a 
D 
3 x 3 
1 4 2 s 40 
2 ? 4 


19 ч ом 2% R 
23K 175 175 325 525 
213 266 29] 277 687 


һ= == 


Ino 


Moyenne 


[2] 
275 


Moyenne 


олз 
0.83 
075 
ЕП 
017 
047 
0.83 
озу 


Moyenne 


675 
E 
2.50 
242 
117 
075 
583 
35 


Écari-T 


601 
1.24 
245 
PAL 
249 
0.90 
БЕ 
146 


Écari-T 


0.66 
1.59 
[IE 
1.00 
ол 
9.38 
204 
0.68 


Écart-T 


6.20 
5.09 
274 
291 
274 
114 
10.12 
37 


p (1) cos wh. 


(8) 


Box et Jenkins (1976) considèrent f (w) comme le spectre de puissance et f* (w) comme 
la fonction de densité spectrale, tandis que KenpaLL (1976) appelle spectre la représentation 


graphique de la densité spectrale f (w) par rapport à о. 
Pour h —0, on obtient par les équations (1)-(2) et (3)-(4) : 


ное: [ ae [ f (0)do 
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et 


»o- | ito | UF (wait = tien | dF (o) =7(0)/с2 = 1 


ce qui montre que la fonction F(w), c'est-à-dire la fonction de distribution spectrale, 
représente bien la décomposition de la variance du processus aléatoire X (г) ct décrit sa 
répartition sur le spectre des fréquences de —z à +n (Jenkins, 1961; CuArFiELD, 1975). 


Les propriétés mathématiques de la fonction F (w) [F* (5] sont celles d'une distribution 
cumulative au sens probabilistique, c'est-à-dire que cette fonction est réelle, monotone, non 
décroissante et bornée, tandis que les propriétés de la fonction f («)[ f * («)] sont celles d'une 
densité de distribution (Jenkins, 1961). Par ailleurs, la fonction de densité spectrale est 
toujours non négative (Davenport, 1970: Агемсотт & Dacuwna-CasreLLE, 1984). 


B) Analyse harmonique 


Tout processus aléatoire stationnaire X (г) auquel correspond une série discrète d'obser- 
vations X,(x,: t— 1, 2, ... №) prises à intervalles réguliers (périodes) peut être représenté 
par un modèle sinusoidal du type : 


X(t) p» cos 01+ sin wi +e(t) (9) 


où p, x et f sont des paramètres estimés à partir des valeurs de la série, w la fréquence 
angulaire et e(t) un processus aléatoire (CuarrigLD, 1975, 1980). 


Le modéle peut étre représenté, dans le cas d'une série finie d'observations, par une 
serie de Fourier (!), c'est-à-dire par une somme de p fonctions élémentaires rigoureusement 
périodiques ( fig. 2) : 


Ni2- 1 
X,—ag- У (a, cos &,1+b, sin 0,1) Һам: COS TI (10) 
E 


oü 
p=1,2 ...N/2-1; ii, РЕ INE 


et о, est la fréquence exprimée en termes d'analyse harmonique ou de la p-ième harmonique, 
c'est-à-dire à, — 2 t p/N et 
а= 
амо 7 E( — ' x/N 
a, 2 (E x, cos о 1)/N 
b, — 2(X x, sin o, 1)/N 


(11) 


Ce modéle est le principe de ce qu'on appelle l'analyse de Fourier ou analyse harmonique, 
les paramètres a, et h, pouvant être estimés par la méthode des moindres carrés. Ils vont 


U} La même sèrie peut étre représentée par un modéle du type autorègressif d'ordre 
+ +т„Х,_„+7„ ой 7, est un processus aléatoire de moyenne u,=0 et de variance с? 
(CuwrriELD, 1975, 1980) 
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х 


fit) ^ 


Nap 


0 100 200 t 
FIG. 2. — La serie des observations X, (A) est décomposée en cinq fonctions élémentaires rigoureusement périodiques 
A=B+C+D+E+F où B—80cosr 1/80, C —40 sin  (/40, D = 20 cos n (/20, E = 10 sin 71/10 et Е = 5cosm (/5. 


définir deux caractéristiques de la p-ieme harmonique : 


l'amplitude, R,= Jaith, (12) 
et la phase, Ф, —tan^ ! ( —b,/a,). (13) 


L'amplitude К, va permettre de définir une nouvelle fonction, 
1(o,) = NR2/(4 n) = N (az -b2)/(4 x) (14) 


dont la représentation graphique par rapport à la fréquence w, est ce que CHaTFIELD (1975, 
1980) appelle le périodogramme; l'aire sous le périodogramme" ainsi défini étant égale à la 
variance totale de la serie d'observations. Par contre, ANDERSON (1971), à la suite de SCHUSTER 
(1898), appelle périodogramme la représentation graphique de l'intensité R? (ou de R, pour 
Schuster) par rapport à la période, T—27/o (cf. Kenbair, 1976) et spectrogramme la 
représentation de R? par rapport à la fréquence, f— 0/(2 x). 

Cette fonction "I(o,). que l'on appellera désormais périodogramme, et que KENDALL 
(1976) appelle intensité du spectre, peut étre estimée directement à partir des valeurs de la 
série d'observations X,(x,:1= 1, 2, N) 


№2 


1 м2 
Lr У x, cos E (> x, Sin @, и) Jo» (15) 


p=i / pz! 


N/2 N/2 
-( X xe) Y xe СЕ 
p) р 
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Снате1еір (1975) montre que l'équation (15) peut s'écrire, 


N-1 


Llo —(co--2 Y, c, cos e, h)/x (16) 
LEX 


où 
N-h 
eum Y (х, 90а /N 
[La 
est l'estimateur de l'autocovariance y (h) et со=52 celui de la variance y(0), ce qui peut 


suggérer que le périodogramme 1 (c,) est un bon estimateur. asymptotiquement non biaisé, 
de la densité spectrale, 


flo, "uoa У y(h) cos oh} (17) 
L h=1 
c'est-à-dire 
lim E[L(o,)] ^ f (w,). (18) 
N =æ 


Toutefois, Í (%,) n'est pas un estimateur cohérent de la densité spectrale / (о). une fois 
que les ordonnées voisines du périodogramme restent intercorrélées, ce qui donne comme 
résultat un périodogramme irrégulier: la non-cohérence doit être éliminée par la suite par 


des procédés de lissage. 


Comme on a vu auparavant la densité spectrale décrit la distribution de la variance du 
processus X (1). Donc. l'aire limitée par la courbe f (w,) et l'axe œ, des fréquences est égale 
à la variance totale c? (fig. 3). Par conséquent. les bandes de fréquence qui correspondent 
aux pics du périodogramme représentent une contribution importante de ces fréquences à 
la variance totale du processus X (1). 


Dans la pratique, la notion de fréquence se réfère à tous les groupes qui peuvent se 
former à partir d'une série X, de N observations, c'est-à-dire aux groupes qui contiennent 
N, N—1. N—2,... 2 observations. Or, les techniques spectrales sont destinées à analyser 
la contribution de chacun de ces groupes d'observations à la variance totale s; de la série X,. 
Leur mérite se comprend mieux si on prend en considération les concepts de l'analyse 
harmonique. En effet, comme il a été signalé précédemment, la série des données X, — qui 
n'est qu'une réalisation du processus X(1) — peut se décomposer en une somme de p 
fonctions élémentaires (harmoniques). 


П a été aussi signalé que le périodogramme 1 (&,) est fonction des paramètres a, et b, 
de ces harmoniques. Par conséquent, les pics du périodogramme dans la fréquence o, 
correspondent aux valeurs élevées des paramètres а, et b. de la 
p-iéme harmonique: autrement dit, leur contribution à la variance de la série X, est impor- 
tante. De plus. tenant compte des harmoniques, la structure de la série devient plus 
compréhensible. En réalité, par l'analyse harmonique, cette structure est décomposée en 
unités (groupes d'observations) qui se répétent périodiquement; autrement dit, il est répondu 
à la question « Après combien d'observations se répéte, en moyenne, le comportement de la 
serie? ». Par la suite, les caractéristiques fondamentales de cette série X, peuvent être 
énumérées et décrites. Par exemple, dans le cas des données du tableau I. un pic possible de 
flop) ou de L(o,) dans la fréquence w,—4 apporte l'information que le comportement de 
la série se répète en moyenne tous les 18 observations (ou tous les 3 mois) et que, par 
consequent. la variation de la série se vérifie principalement dans le domaine temporel. Au 
contraire, un pic dans la fréquence w,= 12 apporte l'information que le comportement de 
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FIG. 3 Spectre de puissance, Ро). ei densité spectrale, /(@,) du processus stationnaire. 
X,205X,-,-025X,-5 -048X,- 37. 


la série se répete toutes les six observations (ou chaque mois) et que la variation se vérifie 
principalement dans le domaine spatial (cf. Annexe I). 


C) Analyse spectrale multidimensionnelle 


L'analyse spectrale peut étre aussi appliquée à la comparaison du comportement de 
deux ou plusieurs séries de données. Ceci est en général plus intéressant que la simple 
analyse d'une seule série de données (PLarr & Denman, 1975). Considérons deux processus 
aléatoires stationnaires, X (1) et Y (1), la fonction d'autocovariance avec retard (LEGENDRE & 
LeGENDRE, 1979) ou autocovariance croisée, Y, (h), est définie par : 

N-h 


y, 0) У IX ()— u][Y (+) VIN, h-0, 1, 2,... N-I (19) 
à laquelle correspond la corrélation avec retard ou corrélation croisée, 


Pas Un = e (Tex (0). Yyy 0) (20) 


Par analogie avec les equations (3) et (7), on définit, en termes de la transformee de 
Fourier, l'autocovariance avec retard ou autocovariance croisée par : 


„®-= | ei** f. (0) do (21) 


-* 
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et la densité co-spectrale par : 
x 
fal] У (De te" (22) 
Km те 

Il faut noter que, tandis que la densité spectrale f (w) était une fonction réelle, la densité 
co-spectrale est une fonction complexe; ceci parce que, contrairement à y (h), y,,(h) n'est 
pas une fonction symétrique (CHatrELD, 1980; KeNpatt, 1976). Ainsi, sachant que 

e '*"'—cos oI—i sin wt, 

fxs (©) est de la forme : 


Fes (090) =k vy (0) — iq, (00) (23) 


dont la partie réelle, le co-spectre, 


k,(0)—[1/(21)) У Yx (h) cos oh (24) 
h 
-[/Q7)] besos Y, [Ye 002 v, 01)] cos zi 
h=1 
mesure la corrélation — en phase — entre les deux processus X (t) et Y (t), et la partie 


imaginaire ou spectre de quadrature, 


409)= [1/02 л)]) У Ya lh) sin oh (25) 


һ= – = 


-[U(Q x)] { У [Ye 00 —,, (А) sin wn} 
h=1 


mesure la corrélation entre les deux séries quand une est décalée par rapport à l’autre de 
1/4 de longueur d'onde (PLarr & Denman, 1975). 


Une autre façon d'exprimer la densité co-spectrale est : 


fey (9) = a, (0) e' * (9 (26) 
où 
ta (0) = kz (®) + qz, (0) (27) 


est l'amplitude du spectre, qui produit des graphiques équivalents à ceux de la densité co- 
spectrale, et 


9,,(0) = tan ! [—4,, (0)/k xy (c)] (28) 


est le spectre de phase de la densité co-spectrale qui mesure la vitesse de variation d'un 
processus par rapport à l'autre (CHATFIELD, 1975, 1980; LEGENDRE & LEGENDRE, 1979). 


Deux autres fonctions sont intéressantes à étudier dans l'analyse spectrale 
multidimensionnelle : la cohérence et la fonction de gain (CuarrieLpD, 1975, 1980; KENDALL, 
1976; LEGENDRE & LEGENDRE, 1979). La cohérence (carré), qui mesure la corrélation linéaire 
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entre les deux processus X (1) et Y (1), est donnée par : 


Ca (0) = [k2, (0) + 42, (0)]/[/,(0)./,(0)) ^ 0SC,, (9) &1 (29) 
=о{, (W) (e). f, (®)] 


ou f, (0) et f, (o) sont les densités spectrales de chacun des processus. La fonction de gain, 
c'est-à-dire le coefficient de régression de Y (t) sur X (г) à la fréquence о, est donnée par : 


Gx (0) — ‚ //,(®). С, (9)/7, (0) (30) 
=% (0)/f, (®). 


En genéral, il est nécessaire que trois de ces paramètres soient pris en considération 
afin que l'analyse fasse apparaitre les interdépendences des deux séries de données. Un trio 
usuel est formé par le co-spectre, la phase et la cohérence. Toutefois, leur interprétation est 
plus difficile que celle de la corrélation croisée (cf. KenpaLz, 1976). Plus précisément, la 
fonction k,,(«) fait apparaitre les liens entre les oscillations harmoniques du premier 
processus avec les oscillations, de méme période, du second. La fonction @,,(w) indique le 
déphasage entre les oscillations de méme fréquence des deux processus, tandis que la fonction 
C,,(«) mesure l'intensité des liens entre les oscillations en question. Une des précautions à 
prendre dans l'estimation du spectre de fréquence est celle de procéder à « l'alignement des 
deux séries de telle facon que la valeur maximale de la fonction de corrélation croisée se 
trouve au voisinage de l'origine » (HuyserecHTs, 1975). 


D) Analyse harmonique multidimensionnelle 
De la méme façon que précédemment [équations (15) et (14)], on peut définir, par 


analogie, le co-périodogramme de deux processus aléatoires X (г) et Y (г) par la fonction 
complexe : 


м2 м2 
өз ( X se) Y пема (31) 
pet pz! 


soit. 
Ixy (Op) — k,,(0,) —1q,, (Op) (32) 
ou 
k sy (0,) = N (apx apy + bpx bp.) /(4 7) (33) 
correspond au co-spectre, 
Axy (Op) = N (ay, b, — ay, b,,)/(4 m) (34) 
au spectre de quadrature, et 
te (0,) = kiy (о) + qz (05) 35) 


a l'amplitude, et où 78 bre 
des processus X (1) et Y (t). 


Comme dans le cas du périodogramme, le co-périodogramme n'est qu'une estimation 
non-cohérente de la densité co-spectrale. 


apy et b, sont les coefficients de Fourier des p-iéme harmoniques 


16 G. P. STAMOU ET J. P. CANCELA DA FONSECA 


Ш. — TECHNIQUES SPECTRALES 
A) Estimation des fonctions spectrales 


Il est évident qu'au fur et à mesure que la valeur de h — décalage ou retard — 
augmente, la précision de l'autocovariance c, diminue. Il semble donc nécessaire de pondérer, 
dans l'estimation de la densité spectrale et du périodogramme correspondant, les valeurs de 
c, par un facteur décroissant À, ou fenêtre avec retard (CHaTFIELD, 1975, 1980; HUYBERECHTS, 
1975) : 


M 


fep (hoc e2 Y, ApC, cos Т, (36) 
h=1 
ou 
M 
fan Y Ac cos у; (37) 
к= см 


où, A, est la fenêtre avec retard, M «NN le point de troncature (« truncation point ») et 
h<M. 


La transformée de Fourier de la fenêtre avec retard est appelée fenétre spectrale 
(Снлаткїкїр‚ 1975, 1980). 


Parmi les fenêtres avec retard, les plus utilisées actuellement sont celles de Tukey et 
de Parzen (Brackman & Tukey, 1959; Parzen, 1961; ANDERSON, 1971; CHATFIELD, 1975; 
Huy8ERECHTS, 1975) z 


(a) Fenêtre de Tukey : 


[1+ со (x h/M)]/2, 0<h<M 
= (38) 
0, h>M 
(b) Fenêtre de Parzen : 
1—6(h/M)?+6(h/M)*,  0<h<M/2 
М= | 2(1—/М)?, M2<h<M (39) 
lo. h>M 


Selon Parr & Denman (1975) la valeur maximale de л doit être choisie entre N/4 et 
N/10. 


CHarreLo (1975, 1980) et Ниувекеснтѕ (1975) ont montré que l'estimateur de la densité 
spectrale fo, ) (éq. 36), comme le périodogramme I(c,) (éqs. 16 et 18), converge vers la 
vraie densité Spactrale quand N => 2 : 


lim E[f («,)] > (9) (40) 


mais qu'il n'est pas, non plus, un estimateur cohérent de cette densité. 


Une autre facon de rendre cohérente l'estimation de la densité spectrale en utilisant 
comme estimateur le périodogramme, T(oy ), est par le procédé de lissage de ce dernier en 
utilisant comme transformation linéaire" ou filtre linéaire les moyennes mobiles (cf. 
Annexe II), c'est-à-dire en groupant leurs valeurs en ensembles de dimension m et en 
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calculant les moyennes respectives : 


Ло) =100,) e (1/m) Y o) (41) 
4 


où, 0,—2nj/N, —m/2<j<m/2 et m est un nombre entier. Le choix de la dimension m 
dépend du type de données et se fait de facon empirique. 


Pour estimer les valeurs de Î(w,) dans l'intervalle [0, (m/2— 1)] ainsi que dans l'intervalle 
[(N/2—m/2 - 1), N/2], la fonction Ѓ(о,) est considérée symétrique en p=0 et р= №2 (voir 
Annexe II). 

Cette méme technique de lissage est utilisée dans l'étude des analyses spectrale et 
harmonique multidimensionnelles, non seulement dans l'estimation des densités spectrales, 
fa (o) et Íy (o), et co-spectrale, p (о), mais aussi dans celle de l'amplitude du spectre, 
$,,(&,), et de sa cohérence, C., (0,) (CHATFIELD, 1975, 1980; LEGENDRE & LEGENDRE, 1979) : 


f. (0) — (1/m) Y. L (oj) (42) 
4 
f, (9) - (Ит) УТ, (o) (43) 
j 
Л (®„)=(1/т) Y. 1, (0) (44) 
1 
&,, (0) = 5, (0,) +43, (a) (45) 
C, (0) 24, (9, f. (0,).f, CA) ^ 0€C, x1 (46) 
ou, 
pim’ 
kolo =N У (a,a,-b, b, )/(4nm) (47) 
1=p-m 
p*m* 
Ge @)=N У (a,b, —a,b,)/(4nm) (48) 
r-p-m* 


etm=2m*+1 


B) Intervalle de confiance de la densité spectrale 


Afin d'évaluer l'amplitude de l'intervalle dans lequel tombent les vraies valeurs de la 
densité spectrale, estimées à partir de N observations, on peut calculer les intervalles de 
confiance de la densité spectrale pour des fréquences différentes (Jenkins & Warrs, 1968: 
CHATFIELD, 1975; Box & Jenkins, 1976). 


. Selon Jenkins et Warrs (1968), tenant compte de l'estimateur de la densité spectrale 
f{w,). la valeur asymptotique 9f (w,)/f (©) suit approximativement la loi du у? avec ё 
degrés de liberté (x2). Le nombre de degré de liberté (de la fenêtre de retard) est donnée 
par. 

M 


8-2N| Y X (49) 


h--M 


Ainsi. la probabilité que la valeur 8f (c,)/f (®„) tombe dans un intervalle de confiance 
de 100 (1— x) % est 


PDá. -a2 €f (®„)// (9) «xi, 43] 2 1—4 (50) 
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les limites de l'intervalle de confiance étant : 
Sf (ә) i-us et M (ogni. us (51) 


Dans trois cas, on n'a pas besoin de calculer le nombre de degrés de liberté par 
l'expression (49). 115 sont donnés directement par les expressions suivantes (CHATFIELD, 1975, 
1980) : 

(a) Fenêtre de Tukey : ӧ= 2,67 N/M. 

(b) Fenêtre de Parzen : 6—3,71 N/M. 

(c) Lissage du périodogramme par groupes de dimension m : 5 =2 т. 


C) Test de signification 


En principe, la simple inspection des graphiques, résultantes de l'analyse spectrale, est 
suffisante pour mettre en évidence les tendances fondamentales de la structure de la série 
des données. D'ailleurs, dans la plupart des cas, les buts de l'étude sont essentiellement 
l'analyse du comportement des séries d'observations. 

Toutefois, parmi les techniques pour tester l'éloignement des conditions d'indépendence 
d'observations également distribuées, le test de BartLeTT (1954) parait assez puissant (RiPLEY. 
1978). П est basé sur le statistique (BARTLETT, 1954; Hannan, 1960) : 


r S м2 


Una] У. oye f о Y. [1(0o)/2x f (o,)] 1 S£N/2 (52) 
pz! 


pat 


Par ailleurs, ce test de signification, non paramétrique, est du type Kolmogorov-Smirnov 
(Hannan, 1960), où 


max JN Upma —$/М/2|<« (53) 


les valeurs de x étant donnees par des tables pour les seuils de signification requis (SIFGEL. 
1956 : tab. M). Les limites de confiance correspondant sont de la forme +4a/ IN, les 
valeurs de ж étant égales à 1,36 et 1.63 pour les seuils de signification (probäbilités) de 
respectivement 0,05 et 0,01 (Hannan, 1960; RieLev, 1978). Ce test, assez précis pour des 
valeurs de N 264, peut être appliqué à n'importe quelle bande de fréquences (RiPLEv, 1978). 

Si, dans un systeme d'axes coordonnées avec S/N/2 en ordonnées et „қ, en abcisses, 
on rapporte les valeurs de Ом; en fonction des valeurs de S/N/2, le graphique resultant 
est une droite passant par l'origine et ayant une pente unitaire (BARTLETT, 1954; HANNAN, 
1960; RipLEv, 1978) 


IV. — DISCUSSION ET APPLICATIONS 


Comme il a été dit dans l'introduction, l'analyse spectrale a été principalement 
développée dans l'étude des séries chronologiques. Néanmoins, elle peut être aussi 
utilisée dans l'étude des séries spatiales, dont le cas particulier des distributions 
(séries) spatiales en treillis. Ici, certaines différences importantes, parfois négligées, 
existent entre les deux types de séries, dont la principale est celle en rapport avec 
le nombre d'observations : dans les séries chronologiques, ce nombre peut tendre 
vers infini dans une seule direction, tandis que dans les séries spatiales en treillis, il 
peut tendre vers infini par différentes voies (TrosrHEIM. 1987). Dans les deux cas, 
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les observations doivent être ordonnées. Ceci pose le problème délicat des séries 
spatio-temporelles. Parce que, si dans les séries chronologiques la succession des 
dates est une relation d'ordre naturelle et dans les séries spatiales une relation 
d'ordre peut être établie entre les différentes observations, dans les séries spatio- 
temporelles, une telle relation est souvent difficile ou impossible à établir (Dunv, 
in lite). 

Il faut cependant faire remarquer que les concepts d'autocovariance et de 
distribution spectrale sont bien définis pour des faibles processus spatiaux station- 
naires (Тзоѕтнегм, 1987). Dans tous les cas, un des principaux problèmes est celui 
de savoir quelle interprétation leur attribuer. Mais, il est évident que l'analyse 
spectrale est d'une grande utilité en suggérant les types de modèle que l'on peut 
construire à partir des données ou ceux qu'il faut éviter (KENDALL, 1976). particuliè- 
rement en ce qui concerne les données obtenues pour des populations naturelles 
(Onp, 1979). 


Bien que trés délicate à manipuler, l'analyse spectrale est néanmoins un outil 
très attrayant et puissant à utiliser dans les études écologiques (LAUREC & BLANC, 
1974; BARTLETT, 1975; PLATT & Denman, 1975; Onp. 1979), où elle présente de 
grands avantages sur l'analyse harmonique « car la quantité importante de bruit 
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FIG 4 — Penodogramme brut des données d' Achipteria holomonensis. A. adultes; B. immatures. 


FIG. 5 -~ Penodogramme lissé (/ = 7) des donnees d'Achipreria holomonensis. A^ adultes. B: immatures 


Tas. HI. = (a) Hlustration du procédé de lissage. Le périodogramme brut (A) est passé successivement par les filtres = 3 (B), m= 5 (C) et 


=7(D) 
A B € ] D 
k 1 = (1, *21,/3 1,+21,+21,)/5 120215 21542177 
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(b) Exemple du calcul de l'estimateur Ї, en fonction des valeurs de 14, . . I; du périodogramme brut. 
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FiG.6 — Périodogramme lissé (т = 3 — 5 — 7) des données d' Achipteria holomonensis. A. adultes; B. immatures. 


fait apparaitre un périodogramme à variance considérable » (IBANEZ, 1972). Dans 
ce type d'analyse, le processus de lissage du périodogramme faisant disparaître des 
détails intéressants (cf. LAUREC & BLANC, 1974), il n'y a pas de règle stricte pour le 
choix des valeurs de m et de M, leurs valeurs étant déterminées par des tentatives 
successives de facon à obtenir un périodogramme qui ne soit pas, ni trop lisse, ni 
trop oscillant. C'est ainsi que CHATFIELD (1980) conseille d'utiliser comme valeurs 
de départ, M égal à 2_ /N et т dans la région de N/40, en particulier dans le cas 
de la fenêtre de Tukey, fenêtre qui a été utilisée dans ce travail. Ainsi, le procédé 
de lissage de la densité spectrale dans l'intervalle (N/20, N/3) va éliminer les pics 
qui apparaissent dans les hautes fréquences, pics qui traduisent la contribution de 
la variation spatiale à la variation totale de la série (CHATFIELD. 1975). П y a donc 
perte d'information spatiale. 


Dans les figures 4 et 5, sont illustrés le périodogramme brut et le périodo- 
gramme lisse (т = 7) des données d'A. holomonensis. Tant le périodogramme brut 
que le périodogramme lissé apparaissent avec un nombre important de pics. donc 
trop inégaux et par conséquent peu explicites. Pour pallier à cet inconvénient, nous 
proposons ici un autre type plus efficace de lissage du périodogramme : le filtrer 
en le passant successivement par trois filtres en série : 


Mag) V (oy) 1” (v, ШЕ] 
—> m3 —> m=5 —> m=] —> 


Ce type de lissage (tab. Ша et IIb; Annexe II), qui n'est qu'une alternative 
des procédés dits de « Hanning » ou de « Hamming » (CHATFIELD, 1975, 1980), est 
une technique beaucoup plus efficace que les autres. Les périodogrammes ainsi 
filtrés apparaissent tres lisses ( fig. 6). 
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FIG. 7 Test de signification de BARTLETT (1954) du périodogramme des données des adultes d' Achiptera holomonen- 
sis. Seuil de signification : p —0.05. 


SIN? 


a NI2 


0 a/ŸN/2 U 
P/N/2 


FIG. 8. — Test de signification de BARTLETT (1954) du periodogramme des données des immatures d'Achipreria holomo- 
nensis. Seuil de signification : p —0.05 
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FIG. 9 Périodogramme des données prefiltrees d' Achipteria holomonensis. A. adultes; B. immatures. 


Les pics dans l'origine, ou plus précisément dans les basses fréquences, tradui- 
sent les apparentes tendances annuelles ou cycles annuels des données (tendances à 
long terme). Notons qu'aprés le procédé de filtrage, les pics ne correspondent 
plus à des fréquences bien définies, mais plutót à des bandes de fréquences. Le 
périodogramme des données qui concerne les adultes présente un maximum 
secondaire dans la bande de fréquences 19-21 ( fig. 6 A). Ce pic, qui correspond à 
un découpage de trois-quatre observations (rhizosphères), reflète la contribution de 
la variation spatiale à la variation totale du recensement. En ce qui concerne les 
immatures, le pic secondaire se situe dans la bande 24-27 ( fig. 6 B), ce qui corres- 
pond à un découpage de trois observations (rhizosphéres). L'allure des périodo- 
grammes — significatifs au seuil de 5% ( fig. 7 et 8) — indique que la variation 
des deux séries de données est principalement temporelle, et secondairement spatiale. 
Les pics secondaires indiquent que l'acarien montre une tendance, bien que faible, 
à former des agrégats plutót autour de la rhizosphére de Festuca sp. qu'autour de 
la rhizosphère d' Arrhenatherum elatius. En effet, la densité (moyenne arithmétique) 
des animaux prélevés autour de la rhizosphére F est supérieure à celle des animaux 
prélevés autour de la rhizosphére A, bien que cette différence ne soit pas significative 
(test-t). 


Il est possible que l'analyse spectrale soit affectée par le fait que les séries 
chronologiques en question soient non stationnaires et présentent une nette tendance 
annuelle. Cette tendance annuelle peut étre, dans les cas moins bien définis, mise 
aussi en évidence par un procédé de lissage, sur un plus ou moins grand intervalle 
de temps (Jacos, 1980), de préférence par les moyennes mobiles (KENDALL, 1976). 
Pour que cette tendance soit éliminée, il est suggéré de filtrer ces séries de données 
avant l'emploi des techniques spectrales (IBANEZ, 1972; LAUREC & BLanc, 1974). 
Nous avons utilisé le filtre le plus simple : 


X'(t) X (t- 1) —X (t), 1=0, 1, ... N-1 (54) 
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FiG. 10 FIG 11 


FIG 10 — Diugrummes des résultats de l'analyse harmonique multidimensionnelle des données des deux series. 


adultes / x) et immatures / y) d' Achipreria holomonensis A, co-spectre: В. spectre de quadrature: C, amplitude du 
spectre: D, spectre de phase; E. cohérence. 

FIG. 11. — Diagrammes des résultats de l'analyse harmonique multidimensionnelle des données préfiltrées des deux 
séries, adultes (x) et immatures (y) d'Achipteria holomonensis. A, co-spectre; B, spectre de quadrature: C, ampli- 
tude du spectre; D, spectre de phase: Е. cohérence 


employé déjà avec succès par LAUREC et BLANC (1974). L'information résiduelle 
apportée par les nouveaux périodogrammes apparait ainsi dans le domaine spatial 
(Ле. 9). 
L'analyse multidimensionnelle des deux séries de données, adultes et imma- 
tures, préfiltrées ou non, confirme ce qui a été constaté précédemment ( fig. 6 et 9). 
L'interdépendance entre les effectifs des adultes et des immatures (co-spectre 
et amplitude; fig. IQA et C) se reflète surtout à l'échelle temporelle (basses 
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FIG. 12 — Periodogramme des données, adultes 4-imma- 
tures, d'Achipteria coleoptratu. A, coucheL: B. 
couche F (ci-dessus) Fu 
FIG. 13 Périodogramme des données d'Achipieria 23 
coleoptrata. А, adultes, couche L; B, immatures. 
couche L: C. adultes, coucheF: D. immatures. - 
couche F (ci-contre) 10 20 30 40 P 


fréquences) et est moindre à l'échelle spatiale (bande de fréquences 20-30). Cette 
interdépendance temporelle est confirmée par les valeurs maximales du spectre de 
quadrature ( fig. 10B et 11 B) et par un déphasage temporel négatif ( fig. 10D et 
11 D). A l'échelle spatiale, l'interdépendance est mise secondairement en évidence 
par le co-spectre et l'amplitude calculés sur les deux séries préfiltrées ( fig. 11 À et 
C). En particulier, le maximum du graphique de phase se situe dans la bande 20- 
30, ce qui indique un déphasage spatial, positif, entre les deux séries, la vitesse de 
variation d'une des séries étant plus importante que celle de l'autre. Pour interpréter 
ce phénomène, une des hypothèses possibles est que les effectifs des adultes oscillent, 
en moyenne, plus rapidement à l'échelle de trois rhizosphéres (F-A-F) que les 
immatures, fait qui refléte la tendance des adultes, par rapport à celle des immatures, 
à former plus d'agrégats autour de la rhizosphére F qu'autour de la rhizosphére A. 


En ce qui concerne la cohérence, les valeurs maximales, accompagnées de 
faibles valeurs de phase, se situe dans la bande de hautes fréquences ( fig. IO E et 
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FIG. 14. Diagrammes des résultats de I lyse harmonique multidimensionnelle des données. adultes +immatures 


d' Achipteria coleoptrata. des deux séries, couche L / x ) et couche F (v). A, co-spectre: В, spectre de quadrature; С. 
amplitude du spectre; D, spectre de phase; E, cohérence. 

FIG. 15. — Diagrammes des résultats de l'analyse harmonique multidimensionnelle des donnees prefiltrees, adultes + im- 
matures d' Achipteria coleoptrata, des deux séries, couche L { x) et couche F / v). A, co-spectre; B, spectre de quadra- 
ture: C, amplitude du spectre: D. spectre de phase: E, cohérence 


11 E). C'est-à-dire, la corrélation linéaire entre les deux séries, adultes et immatures, 
est importante au niveau spatial. Alors, cette corrélation entre les deux séries de 
données s'exprime dans le domaine de deux prélèvements; corrélation qui peut faire 
partie du domaine du hasard. Autrement dit, il existe une relation entre les 
oscillations des adultes et celles des immatures. Il reste à examiner par d'autres 
méthodes, s'il s'agit de relations de type compétition, mutualisme ou autres. 


Si dans le cas d' A. holomonensis on pourrait dire que les techniques spectrales 
apportent une information importante sur la structure des populations, ceci n'est 
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pas le cas de l'emploi de ces techniques aux populations d' A. coleoptrata. Les pics 
des divers graphiques sont assez faibles. Par conséquent, l'organisation des données 
n'est pas tracée avec certitude. Il semble que c'est plutôt le hasard qui se révèle 
comme la caractéristique principale de cette population prélevée, d'ailleurs au 
hasard, à Fontainebleau. Il est possible, en outre, que la faible densité de cette 
espèce affecte gravement les résultats de telles techniques délicates. 


Cependant, méme dans le cas des données d'A. coleoptrata, on peut tirer de 
l'information utile des graphiques spectraux. La variation des effectifs échantillon- 
nés dans la couche superficielle L est tracée principalement à l'échelle de huit 
prélèvements (bande de fréquence voisine de w,=12) (fig. 12). Une variation 
comparable est exprimée par le nombre total d'immatures échantillonnés dans les 
deux couches considérées, L et F (tab. II). Ces observations semblent refléter la 
réalité: par exemple, dans le cas des immatures, le rendement moyen des quatre 
prélèvements faits au milieu de la parcelle est significativement inférieur au rende- 
ment des quatre prélévements faits aux deux bouts (2--2) du biotope (tab. II). De 
plus, les deux graphiques correspondants ( fig. 12 et 13) présentent un maximum 
secondaire dans le domaine spatial de trois prélévements (bande de fréquences 
voisine de Ф, = 32). Or, trois prélèvements ne constituant aucune décomposition 
du nombre total de prélèvements réalisé (8), ceci traduit plutôt l'effet du hasard. 


L'allure des périodogrammes montre que si la variation des effectifs des 
immatures est plutót liée à la couche L, celle des effectifs des adultes l'est avec la 
couche F (fig. 12 et 13). D'une manière générale, les pics des graphiques en 
question se situent dans des fréquences plus ou moins équidistantes, voisines de œ, 
égal à 1, 15, 30 et 45, ce qui reflète plus ou moins l'effet primordial du hasard 
dans la variabilité de ces effectifs. 


L'analyse multidimensionnelle décrivant l'interdépendance des deux séries de 
données, adultes et immatures ( fig. 14 et 15), confirme les conclusions énoncées 
ci-dessus; elle n'apporte donc pas d'information additionnelle. 


V. — CONCLUSION 


Bien qu'on tienne compte des réserves émises sur l'analyse des séries de données 
spatio-temporelles, les mérites des techniques spectrales semblent confirmées dans 
l'analyse des données de recensement d'A. holomonensis, où, d'une part, l'échantil- 
lonnage était défini par des prélévements orientés selon des transects, et, d'autre 
part, la densité de la population était élevée. En effet, par ce type de technique, 
l'information temporelle et l'information spatiale sont analysées en méme temps 
de sorte que seules les caractéristiques fondamentales de la structure des données 
sont révélées. En outre, les paramètres résultant de l'emploi de l'analyse multidimen- 
sionnelle sont tres explicites, notamment ceux de phase et de cohérence. 


Par contre, dans le cas des données de recensement d'A. coleoptrata, où le 
protocole d'échantillonnage avait pour base l'échantillonnage au hasard dans une 
aire déterminée et la densité de la population était faible, l'emploi des techniques 
spectrales n'est pas satisfaisant. Dans ce cas, où les techniques spectrales ne montre 
que l'effet du hasard dans la variabilité des données, les limites de ce type d'analyse 
sont mis en évidence : elles ne doivent pas être appliquées à des données prélevées 
par un échantillonnage au hasard. 
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Toutefois, il faut faire trés attention au type de protocole expérimental choisi, 
en raison du fait, souligné auparavant, que, si la succession des dates est une 
relation d'ordre naturelle, la succession des transects-dates ne l'est pas. Ceci peut 
faire émerger l'influence du temps de facon plus marquante que celle de l'espace 
(Dusv, in litt.). Néanmoins, en biologie du sol, les recensements faunistiques se font 
chaque fois par prélèvement de « carottes » de sol, prélèvement qui va perturber le 
transect échantillonné, celui-ci ne pouvant donc pas être échantillonné la fois 
suivante. L'estimation de l'évolution des populations se fait ainsi à partir d'échantil- 
lons au hasard pris à des endroits successivement différents. Chaque transect étant 
prélevé à une date déterminée, l'ordre temporel est maintenu. Cependant, l'ordre 
spatial, le méme dans chaque transect, bien qu'il puisse paraitre arbitraire, donne 
des séries de valeurs non seulement indicatrices de la variabilité des populations au 
long du transect à un moment donné, mais aussi indicatrices de l'évolution de ces 
populations dans le temps. 


En conclusion, pour que l'analyse spectrale puisse étre appliquée avec toute 
sa puissance, le meilleur échantillonnage doit étre systématique plutót qu'au hasard 
(cf. RiPLEv, 1981) 
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RESUME 


La technique de l'analyse spectrale est décrite avec beaucoup de détail, ainsi qu'une 
nouvelle méthode de lissage des données et des périodogrammes. En attirant l'attention sur 
les difficultes inhérentes à son application aux séries spatio-temporelles, cette technique a 
еге utilisée dans l'étude comparée de la distribution spatio-temporelle de deux espèces 
d'Acariens Oribates du sol, Achipteria holomonensis et Achipteria coleoptrata. La première 
espèce fut obtenue par recensement périodique au long d'un transect dans les rhizosphères 
de deux graminées tandis que la seconde, obtenue aussi par recensement périodique, fut 
échantillonnée au hasard dans la parcelle expérimentale. En conclusion, l'analyse spectrale 
à donné des résultats plus probants pour des données obtenues par échantillonnage systéma- 
tique (données ordonnées) que par un échantillonnage au hasard. 


SUMMARY 


Statistical tools in soil biology 
IX. Spectral analysis and spatio-temporal distribution 


The technique of spectral analysis is described in great detail as well as a new method 
of smoothing the data and the periodograms. Though calling the attention to the inherent 
difficulties of applying this technique to spatio-temporal series, it was applied to the 
comparative study of the spatio-temporal distribution of two species of soil Oribatid Mites. 
The first species has been obtained by periodic census throughout a transect in the rhizosphe- 
res of two grasses whereas the second, obtained also by periodic census, has been sampled 
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at random in the experimental plot. In conclusion, the spectral analysis gave better results 
when the data were obtained by systematic sampling (well-ordered data) rather by random 
sampling 


Annexe I 


Relation entre la fréquence o, et le comportement de la série 


La fréquence o, étant donnée par : w,=2=p/N et la période par Т = 1/0, les valeurs 
de N et de p dans l'exemple du tableau I sont respectivement : 

— N=6 rhizosphères x 12 mois = 72 observations, et 

= р=1, 2, ... N/2=1,2, ... 36. 

Donc, раг conséquence, pour les différentes valeurs de р le comportement de la série se 


répète tous les n observations. Ainsi, pour : 
р=1 w, =0,09 T=11,46 soit 11-12 mois ou 66-72 observations 


p=2 w, 0,7 T= 5,73 soit 6 mois ou 36 observations 

р=3 в, =0,26 = 3,82 soit 4 mois ou 24 observations 

р=4 ©, = 0,35 = 2,86 soit 3 mois ou 18 observations 

p-12  o,,—1,05 = 0,95 soit 1 mois ou 6 observations 

p=24 0@0,,=2,09 T= 0,48 soit <1 mois ou 3 observations 

p=36  @36=3,14 = 0,32 soit <1 mois ou 2 observations 
Annexe II 


Périodogrammes 
Filtres linéaires et lissage par les moyennes mobiles 


La méthode de lissage par les moyennes mobiles donne la tendance de l'évolution des 
séries chronologiques, ainsi que des périodogrammes (voir entre autres, CHATFIELD, 1975, 
1980; KeENpaLL, 1976; LeGENDRE & LEGENDRE, 1979; Orp, 1979). On calcule successivement 
les moyennes arithmétiques de m —2 т* + 1 données contigués que l'on place au centre de 
l'intervalle m des données. 


On transforme ainsi la série initiale, 


Xem X323: Xm [21,2,... N 
dans une autre série, у, par un filtre linéaire, la moyenne mobile, correspondant à 
m —2m* +1 données contiguës 
*m* 
y,-[I/2m*-1) Y x. ой Г=ї+т* 
{=-т* 


Оп peut aussi passer la série initiale par une succession de filtres en calculant successive- 
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ment les moyennes mobiles des m,=2m*+1(r=1,2, ...) valeurs contigués des séries 
obtenues chaque fois (tab. Ша et III b). Donc, par exemple, pour m, —2mf +1. 
.=[1/(2тї+1)] Y ».; où ("tmi 
= m" 


et ainsı successivement. 


De toute façon, pour chaque série de N valeurs, on ne peut calculer que N—2m* 
moyennes mobiles. Toutefois, comme dans le cas du tableau Ша, on a parfois besoin 
d'estimer toutes les valeurs de la nouvelle série correspondant à toutes les valeurs de la série 
initiale. Dans ce cas, on peut attribuer un poids à certaines valeurs de la série initiale pour 
estimer les premiéres et les derniéres valeurs de la nouvelle série qui ne peuvent pas étre 
calculées par la technique habituelle des moyennes mobiles. Pour cela, on considère la série 
symétrique par rapport à ces valeurs. Ainsi, par exemple, dans le tableau Ша, pour 
m=3, Гү=(1-,+1,+1,)/3=(1,+21,)/3, I-> étant égal à L, et pour m=s. 
I; 2( 4T. 10 +13)/5=(1 +21 +215)/5, où T. ; I; et I 5T; 

Exemple du tableau III b : calcul de l'estimateur LU, koe 
x, I. 

Я т=3 I; 2(I5 4-1; -14)/3 

2, m,-3  W-(Ll4I 24 LIES 1/S 

- [IL 19/34- . (1 4-19 -1,9)/3]/5 
=(1,+21,+3 1+3 L 43 1+2 154 1,9)/15 


w, m,=7 =(14+1; +15 +++ ES 401, 9)/7 
=[(1,+2 1,+3 1,43 „+3 1542 L4 L)/I54- ... (L2 144-3 1543 1,543 
1,,21,41,3/15]/7 
=(1,+3 L+6 1,+9 1,+12 154-14 1,+15 1,+14 1,+12 lo+9 T5446 T, +3 
1,5 1,4)/105 
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